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1. Introduction 

Chiral extrapolations of lattice data to the physical pion mass and the continuum or a — > limit 
remain major sources of systematic uncertainty in the determination of hadron masses and matrix 
elements. A test that lattice QCD must successfully pass before predictions can be fully trusted 
is to reproduce known experimental results. One such indicator is the determination of meson 
decay constants, such as f n + and/^+, with phenomenological values of 92.42 ±0.07 ±0.25 MeV 
and 113.0± 1.0 ± 0.3MeV, ^ respectively. The problem is that simulations for smaller quark 
masses rapidly become very costly in computer time. Recent advances have been on two fronts: 
firstly faster machines have become available, with speeds in the Tflop range and secondly the 
hybrid Monte Carlo algorithm used in the simulations has been improved. In particular in the 
new simulations reported here, we have used trajectory length one with three time scales in the 
molecular dynamic step (one for the glue term [Q] and now two |^] for the fermion term in the 
action) which allowed the computationally expensive pieces to be updated less frequently. This 
was coupled with the use of an auxiliary fermion mass, 

The results reported here use Wilson glue (plaquette) and two mass degenerate 0(a) -improved 
Wilson quarks (so effectively we are simulating 2-flavour QCD). As emphasised by Liischer [Q], 
these 'clover' fermions are well understood: in particular the addition of certain irrelevant terms, 
both in the action and operators, and the non-perturbative determination of their coefficients allow 
discretisation errors to be reduced to 0(a 2 ). For example adding the 'clover' term together with 
the appropriate coefficient c sw is sufficient to determine the 0(a) -improved masses, such as the 
pseudoscalar mass m ps while to determine the decay constant, given by 

(0\^\ps) = ^m ps , (1.1) 

the axial current srfn must also be 0(d) -improved, which can be achieved by setting 

^=Z A ^ MP , st™ = (1 + \b A (am qi +am q2 )) (A^ +c A ad^P) , (1.2) 

where A M = q x y M y 5 q 2 and P = q l y 5 q 2 . 

Several years ago we started simulations at four /3 values and reached pseudoscalar masses of 
~ 600 MeV. We have started new simulations at j8 = 5.29 and /3 = 5.40 at lower quark masses. 
Our present status of the lower quark mass runs used in this report is given in table [j] This has 
enabled us to reach pseudoscalar masses of 350 MeV or less. The force-scale was the unit used to 
set the scale, together with a reference value ro = 0.5 fm. Our results for r^/a are shown in fig. [j]. 
In our extrapolations we presently include results for heavier pseudoscalar masses; hopefully the 
situation will improve with more smaller quark mass results being generated so that a linear fit for 
the lighter quark masses will suffice. The extrapolated values of rf^/a in the chiral limit (ro/a) c are 
used to determine the scale. 

2. Chiral perturbation theory 

While the sea quark masses (S) are given (implicitly) in table [T], valence quarks (V) do not 
have to be chosen to have the same mass. Chiral Perturbation Theory, %PT , has been extended to 
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5.29 
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Table 1: Present data sets. The new runs are at (j8, k s ) = (5.29,0.1362), (5.29,0.13632) and (5.40,0.1364), 
where '5" means sea quark. m ps , m^ c are the pseudoscalar and vector particle masses respectively. L is the 
box size. For comparison, experimentally m K + /m p + ~ 0. 18 and m n + ~ 140 MeV. 
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Figure 1: Results for r^/a versus (am s p s s ) 2 for j3 = 5.25 (squares), 5.29 (diamonds) and 5.40 (triangles). 
The new runs are shown in red. Linear fits have been used to extrapolate r^/a to the chiral limit, the results 
being denoted by open symbols. The vertical dashed lines (left to right) represent the chiral limit, and using 
LO jPT approximate positions of a 11, and fictitious ss pseudoscalar particle computed from n + and K + 
respectively. 



Partially Quenched Chiral Perturbation Theory, PQ%PT, [|| 0|. While it is expensive to generate 
dynamical configurations, it is computationally cheaper to evaluate correlation functions on these 
configurations, so that a range of valence quark masses can be used. Using the Leading Order, LO, 
and Next to Leading Order, NLO, results ^ for the pseudoscalar masses and decay constants in 
terms of the quark mass, we eliminate (iteratively) the quark mass from these equations to give for 
degenerate mass valence quarks 

F V P I = fa+h{M S p l) 2 +fc(M V p Jf + f d [{Mf s f + (M^f) H(M s p s s ) 2 + (M v p D 2 ) , (2.1) 

where we have also rescaled the pseudoscalar mass, rrt^f , and decay constant, f£f with say ro c , ie 

M ps = r 0c m ps , F ps = r 0c f ps , (2.2) 
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with A,B £ {V,S}. f a (the LO result) and f, i = b,c,d are given in terms of the low energy 



constants, LECs, ~ —0.76, 0:5 ~ 0.5 (evaluated at a scale /I 
(the decay constant in the chiral limit) [||] by 1 



A x = 4ti f Q ) and f ~ 86.2MeV 



1 



/« — r 0c fo 
f, 

fc 



(47r) 2 r 0c / 
1 



(4^) 2 r 0c / 

f - 1 1 

Jd ~ (47r) 2 r 0f / 4 n/ - 



-n/CS* + -n/ln (2(47rr 0c ./ ) 
1 1 

-a 5 + -?i / ln(2(47rro c /o) 2 ) 



When V = S, eq. (2. 1 ) further simplifies to 



= /« + (A +fc + 2f d ln2)(Mf s f + 2/ d (M^) 2 ln(M^ 



SS\2 
ps ) ' 



(2.3) 



(2.4) 



From eq. (2.1 ) the pion and kaon decay constants can be found. We have two mass degenerate 
sea quarks which we associate with the light quark (I where m/ = (m u + m d ) /2), together with two 
valence quarks, which we associate with either the light, /, quark or the strange, s, quark. Again 
manipulating the structural form of the LO and NLO equations gives the result 

=fa + {fb + fc + 2f d ln2)M 2 + + 2f d M ? - + lnM 2 + 

2 



(2.5) 



fa+ [fb+fd m2 + 




fc+fd\^2- 




-fd 




M% lnM 2 + + ^M 2 n+ In (2M 2 + - M 2 + ) + f d M 2 K+ lnM 2 K+ . (2.6) 



Determining the f a and f, i = b,c,d coefficients means that the pion and kaon decay constants can 
be found. While degenerate quark masses are sufficient, see eq. ( |2. l[ ), for both pion and kaon decay 



constants, only the pion decay constant is possible with just sea quarks, eq. (p.4|). 

Detecting chiral logarithms is a notorious problem, but is necessary as it shows that we are 
entering a regime where #PT is valid. This is particularly difficult for decay constants, as can be 
seen from eq. Q that this term is °c {(M s p s s ) 2 + (M™) 2 ) In ((M s p s s ) 2 + (M™) 2 ) which for fixed 

f VV\2 



[M s p s s ) 2 does not vary much with {M v p J) 2 . 



[^] considering the ratio 



ps '. 

We wish for a term <=c (Af„f ) 2 ln(Afp/) z . As suggested in 
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(2.7) 



-l/(4nf(47tro c fo) 2 ) enhances these chiral logarithms. The disadvantage is that 



(with f d /(n 2 f a ) 

mixed quark mass correlators (V ^ S) must be computed. Note also that eqs. (|2. l[), ( |2.7[ ) probe 
different parts of the #PT expression; as can be seen from eqs. (19) and (20) of [{7|], eq. (2/7) sees 



only the 0(1 /nf) terms, while eq. ( |2.1[ ) probes the remaining 0(1), 0(rif) terms. 



If we had rescaled the pseudoscalar mass and decay constant with r^/a (rather than (ro/a) c as here) r^/a 



(r /a) c (l-r m (M%¥ + . 
fb- 



would just give an additional term — (47r/oro c ) r m in eq. (2.3) in the square brackets for 
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3. Results 



We use the well-established procedure outlined in []|] to compute decay constants. We only 



note here that in eq. (1.2), the improvement coefficient, ca, has been computed non-perturbatively, 
[10], while b& is only known perturbatively (we use a tadpole improved version here, JTT|]). We 
expect, however, that as the quark masses used here are quite small this leads to negligible cor- 
rections. The renormalisation constant has also been non-perturbatively computed, [ |I0[ |TT| ] (the 
differences between these results appear to be 0(a 2 ) and hence vanish in the continuum limit). 

We first investigate to see if we are entering a region where chiral logarithms are becoming 
visible. In fig. |^ we show R defined in eq. (2.7) for j8 =5.29 and K = 0.1359 and 0.1362, together 
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Figure 2: Results for R, eq. (g/7|) for (J8, K s ) = (5.29,0.1359) (circles) and (5.29,0.1362) (squares), against 



(MpJ ) 2 ee (ro c mjJ) . The opaque symbols represent the points where V = S when R = identically. The 
dash-dotted curves are also defined in eq. (2.7) and are plotted with c = —0.01659 (using ro/o ~ 0.218). The 
dashed curve is a fit, yielding c ~ —0.0227 or ro r /o ~ 0.187. Other notation as fig. |l[ 



with the curve also given in eq. (2.7). While we do not expect much influence from the chiral 
logarithm, the curves track the data quite well, indeed out to reasonably large quark masses. So it 
would appear the chiral logarithms are visible of about the expected size. (But note the y-axis scale 
- we have subtracted 1, so really this is a very small effect of 0(1%).) 

To determine the decay constants we must first take the continuum limit of the data, and then 
determine f a and fi, i = b,c,d. But as #PT is an infra-red expansion, while the a 2 — » limit 
is ultra-violet and as we are using 0(a) -improved fermions then we expect that there will be no 
problems with the order of the limits, ie first chiral and then continuum. More drastically we shall 
presently assume that we can ignore any 0(a 2 ) error. This assumption must however be checked 
in the future. 

We now turn to a consideration of the partially quenched results. In fig. || we show the results 
together with a fit from eq. Q2.1[ ). This is a global fit giving one parameter set f a and fj, i = b,c,d 
with values 0.190(16), 0.017(13), 0.030(12), -0.0015(59) respectively. (The fit is reasonably 
good given the fact that the data has three varying parameters: j8, fc s and k v .) We first note the 
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Figure 3: Results for Fj/ = r$ c f£J for the data sets given in table [l]. The fit curves are given by eq. (2.1), 
for a common parameter set. The opaque symbols represent the points where V = S. Other notation as fig^jlj 



the value of f a = ro c /o is in good agreement with the value determined from R. However using 
this value to determine fy (see eq. (2.3)) gives ~ —0.017, indicating that we should be seeing a 
much stronger logarithmic dependence (indeed the curves are almost linear). Furthermore using f a 
and fb, f c gives from eq. (2.3) the values a<\ ~ —0.71, a$ ~ —0.63. 0:4 is in reasonable agreement 
with other phenomenological estimates; but a$ is not. So at the moment there is no unambiguous 
confirmation of #PT. 

Consider now the sea quarks alone. In fig. |] we show these results. The curve joining the 
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Figure 4: Results for F^J = ro c m s p s v for the data sets in table [jj The fit curve (dashed-dotted line) is taken 
from eq. (2^) using the parameters that have been determined from fig. ||. The dashed curves are possible 
phenomenological curves, using eqs. (2.3) and (2.4). The experimental value of F n + = rQf K + is indicated by 
a star and cross for ro = 0.5 fm and 0.467 fm respectively. Other notation as fig. 



points uses the previously determined f a , ft, i = b,c,d coefficients. Consistency is seen. The 
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position of the new results is perhaps surprising because they have dropped to almost below where 
the phenomenological value might lie. Also shown is a possible phenomenological curve from 



eq. ( |2.3| ) for ro = 0.5 fm. As previously found here there is little agreement with the lattice results. 



Reducing the ro scale helps somewhat, the second curve shows the phenomenological results using 



ro = 0.467 fm, a result we estimated previously, see eg (This, of course, has the effect of 

making our pseudoscalar masses larger and box size smaller in table [p.) It would seem that the 
strict applicability of ^PT is restricted to a rather narrow region ro c m s p s s < 1 ; the choice of the 
scale is also rather delicate. Another issue are possible finite-size effects, which we are planning to 
investigate later. 

Finally we note values of f n+ = 77(4) MeV, f K + = 93(1) MeV for r = 0.5 fm and f n + = 
82(5)MeV, f K + = 98(2)MeV for r = 0.467fm. Clearly using f K + to set the scale would, at 
present, make the lattice finer. The dimensionless ratio f K + / f K + is ~ 1.21, 1.19 for ro = 0.5fm, 
0.467 fm respectively, to be compared with {f K + / f n +) e xpt ~ 1.223. 
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